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1.  INTRODUCTION 


DYNACYL  is  a computer  code  used  for  the  solution  of  internal  electro- 
magnetic pulse  (lEMP)  problems  in  right  circular  cylinders;  the  code  includes 
space-charge-limiting  and  air-pressure  effects.  This  class  of  problem  occurs 
typically  when  a source  of  photons  strikes  a cylindrical  body  and  causes 
electrons  to  be  emitted  from  its  surfaces  due  to  photoelectric  and  Compton 
processes.  These  electrons  produced  fields  which,  in  turn,  act  on  the 
charged  particles,  thus  altering  their  trajectories  and  affecting  currents 
and  fields  in  cavities.  When  the  space-charge  fields  are  sufficiently  strong 
to  retard  the  motion  of  electrons  into  the  cavity,  either  by  deflecting 
electrons  radically  to  side  walls  or  turning  them  around  and  propelling  them 
back  to  the  emitted  surface,  then  the  space-charge-limiting  (SCL)  condition 
has  been  realized.  The  electrons  can  also  ionize  gas  molecules  in  cavities 
in  the  system.  When  the  ionization  density  is  high  enough,  migration  of  the 
resulting  light-weight  electrons  away  from  the  heavy  ions  toward  the  cylin- 
der walls  can  reduce  the  space-charge  barrier  fields  and  permit  greater 
penetration  of  primary  electrons  into  the  space,  thereby,  significantly 
changing  currents  and  fields  in  the  cavities.  The  goal  of  an  lEMP  analysis 
is  to  determine  currents  reaching  vital  components.  In  this  regard,  both 
space-charge-limiting  and  air-pressure  effects  may  be  important.  DYNACYL 
is  designed  to  aid  in  these  analyses  by  treating  both  effects  simultaneously. 

This  report  contains  a broad  overview  of  the  code;  details  on  how  to 
use  it  are  emphasized.  Detailed  descriptions  of  the  physics  and  modeling 
employed  in  the  code  are  found  in  the  references.  Here,  specific  input 
card  Information  and  output  descriptions  are  given.  Computer  requirements 
are  outlined,  a sample  calculation  Is  discussed,  and  other  practical 
Information  necessary  for  obtaining  useful  Information  from  DYNACYL  is 
given. 


DYNACYL  contains  as  a subset  its  predecessor  code  TEDIEM-RZ.  This 
earlier  version  is  a quasi-static  solution  of  the  same  configuration  as 
that  treated  by  the  dynamic  DYNACYL  segment,  but  does  not  treat  air 
pressure  effects.  It  can  be  invoked  following  directions  of  this  manual, 
but  is  otherwise  left  out  of  discussions.  Complete  documentation  and  a 
user's  manual  for  this  older  code  are  available  in  Reference  1. 


2.  DYNACYL  DESCRIPTION 

DYNACYL  solves  the  two-dimensional,  time-dependent  lEMP  problem  for 
a cylinder  whose  end-plates  emit  electrons  into  the  interior  in  axisymmetric 
distributions.  Particles  representing  large  numbers  of  emitted  electrons 
are  formed  and  injected  into  the  spatial  grid  at  various  energies,  angles, 
and  positions,  and  with  various  amounts  of  charge  depending  on  the  pulse 
height  at  the  time  of  emission.  All  of  the  electron  emission  information 
must  be  specified  to  the  code.  DYNACYL  automatically  provides  for  the  de- 
lay of  electron  emission  from  the  back  face  by  an  amount  equal  to  the  flight 
time  across  the  cavity  length. 

The  particles  of  charge  are  moved  via  Newtonian  equations  of  motion. 
They  are  converted  to  current  densities  on  the  spatial  cell  boundaries  and 
also  ionize  air  molecules  if  such  are  present  in  the  cylinder.  The  current 
densities  are  used  as  source  terms  for  Maxwell's  equations,  and  the  ioniza- 
tion rates  are  used  in  the  plasma  solution  for  the  background  gas.  Both 
Maxwell's  equations  and  the  plasma  condition  are  solved  using  finite- 
differencing  techniques.  The  plasma  treatment  couples  into  the  field  equa- 
tions by  means  of  effective  conductivity  in  each  zone  due  to  the  secondary- 
electron  mobility. 

The  details  of  the  physics  and  modeling  can  be  found  in  Reference  2. 

A useful  discussion  of  the  range  of  validity  of  the  plasma  formulation 
and  other  DYNACYL  treatments  for  a specific  lEMP  air-pressure-effects 
problem  are  found  in  Reference  3.  A comparison  of  code  performance  in 
dynamic  versus  quasi-static  lEMP  problems  is  found  in  Reference  4.  It 

2.  T.  N.  Delmer  et  al.,  "SGEMP  Phenomenology  and  Computer  Code 
Development,"  DNA  3653F,  November  11,  1974. 

3.  1).  C.  Osborne,  R.  H.  Stahl,  and  T.  N.  Delmer,  "Large-Area  Electron- 
Beam  Experiments,"  INTEL-RT  8101-011,  July  15,  1975. 

4.  E.  P.  Wenaas  and  A.  J.  Woods,  "Comparisons  of  Quasi-Static  and 
Fully  Dynamic  Solutions  for  Electromagnetic  Field  Calculations  in  a 
Cylindrical  Cavity,"  IEEE  Trans.  Nucl.  Sci.  NS-21,  December  1974. 
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should  be  noted  that  DYNACYL  has  no  provision  for  field  acceleration  of 
secondary  electrons;  hence  the  code  cannot  treat  the  occurrence  of  avalanch 
ionization. 

A flow  chart  summarizing  the  DYNACYL  code  cycle  is  shown  in  Figure  1. 
Notice  that  electron  emission  information  comes  from  a source  external  to 
DYNACYL. 
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Figure  1.  DYNACYL  flow  chart 


3.  DYNACYL  COMPUTER  REQUIREMENTS 

DYNACYL  is  a FORTRAN-IV  computer  program  of  about  5500  cards, 
including  both  the  dynamic  and  quasi-static  segments  of  the  code  and 
peripheral  plot  routines.  No  machine  language  coding  is  employed.  The 
current  version  operates  on  the  CDC  7600  computer,  although  an  older 
UNIVAC  1108  version  exists.  The  older  edition  is  not  considered  here. 
Core  requirements  are  65000j^q  words  of  small  core  and  100,000^^  words 
of  large  core.  Three  fast-access  files  are  also  required  during 
execution,  one  for  plot  storage  and  two  for  particle  information  storage. 
The  RUN  compiler  is  currently  used,  although  conversion  to  the  faster 
FTN  compiler  should  be  straightforward  at  this  time.  Run  times  vary 
from  1 to  30  minutes  central  processor  time,  depending  on  problem 
conditions.  Core  must  be  preset  to  zero. 
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4.  DESCRIPTION  OF  INPUTS 


4.1  INTRODUCTION 

Detailed  descriptions  of  DYNACYL  input  cards  are  given  in  this 
section.  The  goal  here  is  to  explain  how  to  specify  the  various  physical 
and  computational  parameters  to  the  code.  A brief  definition  section  is 
given  in  Table  1,  in  which  certain  terms  peculiar  to  the  DYNACYL  user's 
manual  and  code  outputs  are  defined.  These  descriptions  are  not  of 
specific  code  variables,  but  rather  of  general  terms  employed  in  modeling 
and  interpreting  problems.  Specific  variable  descriptions  are  found  in 
the  Input  Cards  and  Sample  Problem  sections  below. 

Detailed  discussions  of  card  inputs  are  given  including  variable  name 
and  definition,  minimum  and  maximum  allowed  values  and  numbers  of  values, 
defaults,  locations  on  the  cards,  etc.  Numerical  limitations  on  code 
input  and  calculational  variables  are  summarized  in  Table  2. 

4 . 2 INPUT  CARDS 

Input  card  details  arc  given  in  Table  5.  Several  conventions  are  used  to 
help  in  quick  scanning  of  the  inputs.  If  the  word  "DEBUG"  is  found  at  the 
beginning  of  a variable  description,  that  variable  is  used  generally  for 
debug  purposes  only  and  may  be  ignored  for  standard  code  calculations. 

The  word  "EDIT"  indicates  that  the  variable  is  used  in  choosing  desired 
code  outputs.  Most  of  these  variable  types  appear  on  the  option  card 
(card  2) . Variables  without  the  "EDIT"  or  "DEBUG"  notation  are  generally 
the  physical  inputs  to  the  problem. 

The  designation  "card  number"  actually  means  "card  type".  If  an 
array  of  values  is  bej.ng  read  in  requiring  more  than  one  card,  the  cards 
have  the  same  card  number  until  the  desired  number  of  values  are  specified. 
Those  input  variables  which  are  arrays  and  may  require  multiple  values  can 
be  determined  quickly  by  noting  the  format  specification,  where  multiple 
values  are  indicated  by  a number  appearing  before  the  alphabetic  character 
giving  the  variable  type.  Also,  the  descriptions  tell  how  many  variables 
to  read  in. 
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Table  1 


j 

DEFINITIONS  OF  TERMS  FOUND  IN  THIS  USER'S 
MANUAL  AND  IN  DYNACYL  OUTPUT 


Forward  emission  face: 

Cylinder  face  from  which  forward 
emission  of  electrons  occurs.  Usually 
taken  to  be  bottom  (z=0)  face  of 
cylinder.  Emission  always  occurs 
from  this  face  in  the  code.  If  back 
emission  is  done,  opposite  face  emits, 
also . 

Back  of  cylinder: 

Face  away  from  forward  emission  face 
(z=L) 

L: 

Cylinder  length  (cm) 

Top  of  cylinder: 

Face  away  from  forward  emission  face 
(z=L) 

Bottom  of  cylinder: 

Forward  emission  face  (z=0) 

R: 

Cylinder  radius  (cm) 

D: 

Cylinder  diameter  (cm) 

Mini-print : 

Short  printout  (6  lines)  containing 
summary  and  plot  information  useful 
in  interpreting  calculational  quality 
and  code  output.  Can  be  specified 
at  times  distinct  from  the  large 

2-D  prints . 

2-D  print: 

Large  printout  containing  spatial 
distributions  of  most  physical 
quantities . 

Primary  electrons: 

Designates  electrons  emitted  into  the 
cylinder  by  external  source  such  as 
photo-electric  emission. 

Secondary  electrons: 

Designates  electrons  caused  by 
ionization  of  the  background  gas  by 
primary  electrons.  Modeled  as  a 
plasma  in  DYNACYL. 
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Table  2 


MINIMUM  AND  MAXIMUM  ALLOWED  VALUES 
OF  DYNACYL  INPUT  AND  CALCULATIONAL  VARIABLES 


Defaults  are  given  where  applicable. 


Quantity 

Minimum 

Value 

Allowed 

Maximum 

Value 

Allowed  Default 

Radial  zones : 

(NPG) 

3 

20 

Axial  zones : 

(NRGP) 

3 

50 

Emission  Specifications 

Angular  bins: 

1 

20 

Polar: 

(NANG) 

1 

20 

Azimuthal: 

(KOPT(IO)) 

1 

No  limit  except  by  number  of  particlei 
allowed  to  be  emitted  per  time  step 
(below) 

Energy  bins : 

1 

25  (15  for  time-dependent  spectra) 

Positions : 

1 

20  Number  of  radial  zones 

No.  of  different  spectra: 
(KOPT(38)) 

1 

10  1 

Number  of  different 

particle  time  steps: 
(K0PT(12)) 

Number  of  point  pairs 
defining  emission  pulse 

1 

10 

shape : 

(K0PT(13)) 

2 

41 

Number  of  particles 

emitted  per  time  step: 

Number  of  particles 
being  followed  in  a given 

0 

2000 

time  step: 

Number  of  point  pairs  on 
field  & current  time 

0 

No  limit  except  computer  disk  space 

history  plots: 

2 

1200 

1 1 


Tabic  3. 
INPUT  CARDS 


Card  Columns 

Number  (Format)  Variable  Description 


1 

1-78 

(13A6) 

TITLE 

Title 

2 

1-2 

KOPT(l) 

EDIT: 

(12) 

(I PLOT) 

0 

1 

2 

2 

3-4 

K0PT(2) 

EDIT: 

(12) 

2 

5-6 

K0PT(3) 

EDIT 

(12) 

2 

7-8 

K0PT(4) 

(12) 

(NESKIP) 

2 

9-10 

K0PT(5) 

EDIT 

(12) 

2 

11-12 

(12) 

K0PT(6) 

DEBUG 

2 

13-14 

(12) 

K0PT(7) 

EDIT: 

Card 


Particle  trajectory  information 
No  information 

Plot  sample  trajectories  for  time 
of  DTF  nsec  (see  card  11) 

Print  velocities  and  fields  acting 
on  particles  over  the  trajectories 
as  well  as  plotting  them 

Print  out  every  K0PT(2)  values  of 
the  array  of  emission  particle  data 
(initial  position,  velocity,  charge, 
etc.)  if  zero,  code  calculates  value 
such  that  about  50  sample  particles 
are  printed 

If  >0,  print  out  emission  charge 
characteristics  as  in  K0PT(2)  only 
for  2 X KOPT(38)  different  times. 

Useful  in  time-dependent  spectrum 
calculations  for  verification  of 
emission  data  at  different  times. 

Input  emission  electrons  every  K0PT(4) 
particle  time  steps. 

If  >0,  calculate  averages  and  standard 
deviations  of  certain  summary  quant ies 
(see  mini-print  description)  every  K0PT(5) 
particle  time  steps 

Particle  path  information  for  sample 
trajectories  which  is  printed  out  as 
the  calculation  unfolds. 

Print  fields  and  currents  for  every 
spatial  zone  away  from  forward  emission 
face  up  to  and  including  the  K0PT(7)th 
zone.  See  KOPT(20).  DEFAULT  - 1 


, ) 


Card 

Number 


2 


2 


2 


2 


2 


2 


2 


2 


Table  a.  (cont.) 

Columns 

(Format)  Variable  Description 


15-16 

(12) 


17-18 

(12) 


19-20 

(12) 


21-22 

(12) 


23-24 

(12) 


25-26 

(12) 


27-28 

(12) 


29-30 

(12) 


KOPT(8)  0 Emission  current  density  uniform  radially 

across  emission  face  (s) 

1 Emission  current  density  falls  off  radi- 

ally to  a value  of  1/BFALL  (see  card  11), 
its  value  on  the  cylinder  axis. 


KOPT(9)  DEBUG:  Print  out  information  on  2-D  prints  in 

addition  to  standard  fields  and  currents. 
(B,  charge  per  zone,  potential  relative 
to  side  wall,  etc.) 

KOPT(IO)  Number  of  emission  electron  azimuthal 

(NPHI)  angular  bins 

MIN  = 1 
MAX  = 2000 


KOPT(ll)  Number  of  emission  electron  energy 

(NSPD)  bins.  Cards  13  and  14  read  in  if 

KOPT(ll)  >0.  Values  of  -1  and  <-l 
can  be  used  where  code  sets  mono- 
energetic  or  analytic  spectra  for 
debug  purposes  (see  subroutine  TDIMR7, ) 
MIN  NO  = 1 

MAX  NO  = 25  (15  for  time  dependent 
spectra;  see  KOPT(38)) 


K0PT(12)  Number  of  different  particle  time  steps 

(NTSTEP)  (read  in  on  cards  7 and  8) 

MIN  = 1 

MAX  = 10 

K0PT(13)  Number  of  point  pairs  used  to  define 

(NTPULS)  arbitrary  emission  current  pulse 

shape  read  in  on  cards  9 and  10. 

MIN  NO  » 2 
MAX  NO  - 41 


KOPT(14)  Number  of  radial  positions  at  which 

(NEPTS)  charge  particles  are  emitted 

DEFAULT  = No.  of  RADIAL  ZONES 

MIN  = 1 

MAX  - 20 

K0PT(15)  EDIT:  Print  2-D  prints  every  K0PT(15)  particle 

(NTSKIP)  time  steps.  Overridden  by  DTPRNT 

(Card  6)  if  it  is  non-zero. 


13 


I 


i 


Card 

Number 

Columns 

(Format) 

Table 

Variable 

■3. 

[cont . ) 

Description 

2 

31-32 

(12) 

K0PT(16) 

DEBUG 

: Print  2-D  prints  every  particle  time  step 
up  to  and  including  step  no.  K0PT(16)  if 
>0. 

2 

33-34 

(12) 

K0PT(17) 

EDIT: 

Number  of  different  plot  files  to 
overlay  from  results  of  this  computer  run. 
An  unlimited  number  of  input  decks  can  be 
stacked,  but  only  up  to  the  first  3 can 
be  overlayed 

MIN  = 0 

MAX  = 3 

2 

35-36 

(12) 

K0PT(18) 

EDIT: 

If  >0,  all  2-D  prints  suppressed  until 
more  than  K0PT(18)  particle  steps  have 
elapsed . 

2 

35-36 

(12) 

K0PT(19) 

DEBUG 

: If  >0,  positions,  velocities,  etc.  of 
every  particle  are  printed  out  every 
time  step.  Uses  much  paper. 

2 

39-40 

(12) 

KOPT(20) 

EDIT 

: Print  fields  & currents  at  every  K0PT(20) 
spatial  grid  positions  in  2-D  prints. 
Default  gives  approximately  10  zones  in 
each  direction.  See  K0PT(7) 

2 

41-42 

(12) 

K0PT(21) 

DEBUG 

: Provides  additional  information  on  2-D 
prints  if  >0 

2 

43-44 

(12) 

K0PT(22) 

EDIT: 

If  >0,  plot  total  currents  striking 
bottom,  side,  and  top  on  same  graph  with 
emission  current  wave  form.  All  curves 
are  divided  by  peak  emission  current. 

2 

45-46 

(12) 

KOPT(23) 

EDIT: 

If  >0,  plot  peak  potential  on  cylinder 
axis  versus  time.  The  position  of  the 
peak  can  vary  with  time. 

2 

47-48 

(12) 

K0PT(24) 

EDIT: 

If  >0,  plot  average  and  peak  emission 
electron  energies  versus  time  if  the 
spectrum  is  time-varying  (see  K0PT(38)) 

2 

49-50 

(12) 

K0PT(25) 

EDIT: 

If  >0,  plot  total  current  collected  by 
top  of  cylinder  (opposite  forward  emission 
face)  versus  time  both  before  and  after 
smoothing  numerical  noise.  Plots  every 
time  step  up  to  number  1200.  Current 
is  divided  by  peak  emission  current. 
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Table  5.  (cont . ) 


i 


Card 

Number 

Columns 

(Format) 

Variable 

Description 

2 

51-52 

(12) 

KOPT(26) 

EDIT: 

Same  as  KOPT(25)  only  for  side  of 
cylinder . 

2 

53-54 

(12) 

KOPT(27) 

EDIT: 

If  >0,  plot  total  current  collected  by 
top  of  cylinder  (opposite  forward 
emission  face)  at  radius  less  than  Ro 
and  at  radius  greater  than  Ro  where 

Ro  = FCOLT  * R2.  R2  is  cylinder 
radius  and  FCOLT  is  set  at  beginning  of 
subroutine  TDIMRZ  (currently  set  to  .67). 

Current  is  divided  by  peak  emission 
current. 

! 

2 

55-56 

(12) 

KOPT(28) 

EDIT: 

Same  as  KOPT(27)  only  for  bottom  of 
cylinder  (forward  emission  face). 

i 

\ 

j 

2 

57-58 

(12) 

KOPT(29) 

EDIT: 

If  >0,  plot  axial  electric  field  at  i 

forward  emission  face  (z  = 0)  on  axis 

and  at  half  cylinder  radius  (r/2)  as 

well  as  peak  radial  electric  field 

anywhere  along  cylinder  side  versus  time. 

Position  of  radial  field  can  vary  i 

with  time.  i 

2 

59-60 

(12) 

KOPT(30) 

EDIT: 

If  >0,  plot  axial  electric  field  on 
axis  at  face  opposite  forward  emission 

face  (z  ■ L)  versus  time  1 

• 

2 

61-62 

(12) 

K0PT(31) 

EDIT: 

j 

If  >0,  plot  magnetic  field  at  back  face  j 

(z  “ L)  at  side  of  cylinder  (r  ■ R) 
versus  time. 

2 

63-64 

(12) 

KOPT(32) 

EDIT: 

Input  check-plot  emission  current 
pulse  shape  if  >0. 

2 

■ 

65-66 

(12) 

KOPT(33) 

EDIT. 

If  >0,  plot  spatial  distribution  of  i 
electric  potential  on  cylinder  axis  1 
at  each  2-D  print  time  | 

: 1 

' 

1 

67-68 

(12) 

KOPT(34) 

EDIT: 

1 

If  >0,  plot  net  charge  density  (ions  1 

and  electrons)  spatial  distributions  ! 

for  several  radial  positions  each  2-D  i 

print  time  | 

[ 

1 

} 

t 

1 

1 

I 
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Card 

Number 

Columns 

(Format) 

Table  3. 

Variable 

(cont . ) 

Description 

2 

60-70 

(12) 

KOPT(35) 

EDIT: 

Print  mini-print  every  K0PT(35)  particle 
time  steps. 

DEFAULT  = No  mini-print 

2 

71-72 

(12) 

KOPT(36) 

EDIT: 

Print  pressure  effects  quantities  in 

2-D  prints  in  addition  to  fields  and 
currents.  Includes  conductivities, 
conduction  currents,  charge  densities, 
ionization  rates 

2 

73-74 

(12) 

KOPT(37) 

DEBUG 

: Allows  different  terms  in  particle 
equations  of  motion  to  be  tested  in- 
dependently: 

0 

E and  B both  act  on  particles 

1 

E only  acts  on  particles 

2 

B only  acts  on  particles 

3 

No  force  due  to  fields  acts  on  particles 

In  all  of  the  options  above,  the  slowing 
down  of  particles  due  to  ionization  of 
the  gas  by  the  primaries  still  occurs. 

2 

75-76 

(12) 

KOPT(38) 

Number  of  different  spectra  repre- 
senting complete  emission  current  pulse. 

If  more  than  1 spectrum  is  used,  code 
linearly  interpolates  in  time  from  one 
spectrum  to  the  next.  If  0,  random- 
izing routines  (problem  specific)  are 
called  which  operate  on  the  input  dis- 
tributions to  randomize  the  energy  and 
angular  spectra  (see  subroutines  RAMSPD 
and  RAMIT) . In  that  case,  the  number 
of  spectra  is  the  absolute  value  of 
KOPT(38) 

DEFAULT  "=  1 
MAX  - 10 


2 

77-78 

(12) 

KOPT(39) 

Quasi-Static  only:  read  Green's  function 
for  fields  from  input  unit  KOPT(39) 

2 

79-80 

(12) 

K0PT(40) 

Quasi-Static  only: 
for  fields  to  unit 

write  Green's  function 
K0PT(40) 
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Card 

Number 

Columns 

(Format) 

Table  3. 

Variable 

(cont . ) 

Description 

3 

17-25 

(E9) 

PATH 

Gas  pressure  in  cylinder  in  atm.  Also 
used  to  flag  quasi-static  calculations 
which  can  be  done  for  vacuum  only.  Set 
to  -l.E-20  to  indicate  quasi-static. 

3 

26-34 

(E9) 

R2 

Cylinder  radius  (cm) 

3 

44-52 

Z2 

Cylinder  length  (cm) 

3 

53-61 

(E9) 

NRGP 

No.  of  spatial  zones  in  radial  direction 

3 

62-70 

NPG 

Number  of  spatial  zones  in  axial  direction 

3 

71-79 

(E9) 

NANG 

Number  of  polar  angular  bins  for  emission 
current.  If  >0,  code  reads  in  angular 
distribution  on  cards  15  and  16  with 

NANG  bins.  If  <0,  code  set  cosine  9 
(measured  from  surface  normal)  angular 
distribution  with  NANG  bins  and  cards 

15  and  16  are  not  read  in.  If  <-100, 
code  sets  isotropic  angular  distribution 
with  NANG  + 100  bins  and  cards  14  and 

15  are  not  read  in. 

4 

1-40 

(ID 

LRATIO 

Determines  radial  zoning  in  quasi-static 
case.  Read  in  one  LRATIO  value  for  each 
radial  zone.  The  cylinder  radius  is  sub- 
divided into  minimum  increments  equal  to 
the  radius  divided  by  the  sum  of  the  LRATIO  * 

values.  The  radial  zones  are  then  the 
product  of  the  individual  LRATIOs  times 
the  minimum  increment.  Permits  rapid 
particle  finding.  Read  in  as  blank  in 
dynamic  calculations  (PATH  ^ 0) . 

5 

1-40 

(11) 

LZRAT 

Same  as  card  4 only  for  axial  zones 

6 

17-25 

(E9) 

TMAX 

Maximum  simulation  time  (nsec) 

i 


Card 

Number 

Columns 

(Format) 

Table  3. 

Variable 

(cont . ) 

Description 

6 

26-34 

(E9) 

BAKEMT 

Back  emission  fraction.  If  >0,  code 
emits  particles  from  back  face  of 
cylinder  (Z  = L)  with  same  spectrum 
and  pulse  shape  as  forward  emission. 

Peak  emission  current  density  equals 

BAKEMT  times  peak  forward  emission 
current  density  (card  11).  Back 
emission  is  delayed  by  photon  flight 
time  to  back  face. 

6 

35-43 

(E9) 

SIDEMT 

Fractional  emission  from  side  of 
cylinder.  Not  operational. 

6 

44-52 

(E9) 

REFQ 

Charge  reflection  fraction  for  elec- 
tron back-scattering  simulation.  Not 
checked  out  thouroughly.  If  0,  no 
backscatter ing  occurs. 

6 

53-61 

(E9) 

FINJEK 

Emission  current  radial  extent  is  from 

0 to  FINJEK  X cylinder  radius. 

DEFAULT  = 1.  See  BFALL  (card  11) 

6 

62-70 

(E9) 

DTPRNT 

Print  out  2-D  prints  every  DTPRNT 
nsec  if  >0.  Overrides  K0PT(15) . 

7 

17-70 

(7E9) 

TSTEP 

Time  range  edges  in  nsec  over  which 
particle  time  step  is  constant  and 

defined  by  DTSTEP  (card  8).  ie: 
particle  time 

defined  by  DTSTEP  (card  8),  ie, 
particle  time  step  = DTSTEP  (I)  for 
TSTEP  (I)  < TIME  < TSTEP  (I  + 1) . 

Read  in  KO^  (12)  values.  First  one 
must  be  zero.  TSTEP  (KOPT  (12)  + 1) 
set  to 

MIN  NO  = 0 
MAX  NO  = 10 

8 17-70  DTSTEP  Particle  time  steps  in  nsec.  See  card 

(7E9)  7 above.  If  DTSTEP  (1)  = 0.,  code  sets 

particle  time  step  so  an  average-energy 
particle  traverses  one  axial  zone  per 
particle  time  step.  If  DTSTEP(l)  = -1., 
code  makes  sure  time  step  is  at  least 
less  than  or  equal  to  1/5  emission 
current  pulse  rise  time. 
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Card 

Number 

Columns 

(Format) 

Tal)le  5. 

Variable 

( font . ) 

Description 

9 

17-70 

(7E9) 

TPULSE 

Times  at  which  emission  current  pulse 
defined  (see  card  10).  Code  linearly 
interpolates  between  values  for  emission 
current  level.  Read  in  K0PT(13)  values. 
MIN  NO  = 2 

MAX  NO  = 41 

10 

17-70 

(7E9) 

PULSE 

Relative  emission  current  pulse  heights 
at  times  defined  by  TPULSE  (card  9). 

Read  in  K0PT(13)  values.  Code  shuts 
off  pulse  immediately  after  last  value 
of  TPULSE. 

MIN  NO  = 2 

MAX  NO  = 41 

11 

17-25 

(E9) 

NTOT 

2 

If  >0.,  number  of  electrons/cm  emitted 
at  axis  over  entire  emission  current 
pulse.  See  BFALL  if  radial  variation  of 
emission  current  desired.  Either  NTOT 
or  JPEAK  must  be  Input. 

11 

26-34 

(E9) 

JPEAK 

If  >0.,  peak  emission  current  density 
in  amp/cm^  at  axis  of  forward  face.  If 
zero,  calculated  from  NTOT  above. 

Either  JPEAK  or  NTOT  must  be  Input. 

See  BFALL  If  radial  fall-off  of  emission 
current  desired. 

11  35-43  DTF  EDIT:  Time  Interval  In  nsec  over  which  sample 

(E9)  particle  trajectories  are  plotted  (see 

KOPT(l)).  Particles  are  plotted  for 
DTF  nsec,  then  a new  batch  is  tagged 
and  followed  for  the  next  DTF  nsec 
throughout  the  entire  problem. 

DEFAULT  = average  energy  particle  time 

to  transit  the  cavity  assuming 
no  SCL 

11  44-52  BFALL  Emission  current  radial  fall-off  if 

(E9)  K0PT(8)  set  to  1.  Emission  current  falls 

off  to  a value  at  cylinder  outer  wall  of 
1. /BFALL  times  its  value  on-axis.  See 
also  FINJEK  (card  6).  This  option  de- 
termines radial  fall-off  rate,  but  FINJEK 
determines  extent  of  emission  region 
Independently. 


Card 

Number 


12 


13 


14 


15 


16 


Table  3.  (cent.) 

Columns 

(Format)  Variable  Description 


17-70  TSPECT 

(7E9) 


1-80  EEDGE 

(8E10) 


Read  in  only  if  K0PT(38)  >1  for  time- 
dependent  emission  current  spectra.  These 
are  the  times  in  nsec  corresponding  to  the 
input  spectra  on  cards  13  and  14.  Each 
set  of  cards  13  and  14  corresponds  to 
One  value  of  TSPECT.  For  mono-energetic, 
time-varying  emission  currents,  replace 
cards  13  and  14  with  a card(s)  (16X,  7E9) 
On  which  emission  electron  energies  in 
kev  are  specified  corresponding  to  the 
times  on  card  12. 

DEFAULT  = 0 
MAX  NO  = 10 

Emission  electron  spectrum  energy  bin 
edges  in  kev.  See  card  14.  Read  in 
KOPT(ll)  + 1 values  if  KOPT  (11)  > 0. 


1-80  HELN 

(8E10) 


1-80  ANGEDG 

(8E10) 


1-80  NTHETA 

(8E10) 


Emission  electron  spectrum  relative  no. 
of  electrons  per  unit  energy  in  energy 
bins  defined  by  card  13.  Read  in  KOPT(ll) 
values  if  KOPT  (11)  > 0.  For  time- 
dependent  spectra,  KOPT  (38)  pairs  of 
cards  13  and  14  are  read  in. 

Emission  electron  angular  distribution 
bin  edges  in  degrees.  Measured  from 
Surface  normal.  Read  in  NANG  + 1 values 
(see  card  3)  if  NANG  >0,  otherwise  omit. 

Emission  electron  angular  distribution 
relative  number  of  electrons  per  unit 
solid  angle  in  bins  defined  by  card 
15.  Read  in  NANG  value  if  NANG  > 0, 
otherwise  omit  (see  card  3)  . 
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5.  SAMPLE  PROBLEM 

A sample  problem  is  described  in  this  section  in  which  code  inputs 

and  selected  outputs  are  explained.  The  problem  modeled  is  a cylinder 

-4 

15  cm  long  by  30  cm  diameter  containing  air  at  pressure  2.63  x 10  atm 
(200  pm  Hg) . A pulse  of  electrons  is  emitted  perpendicularly  and  uniformly 
fromoneend  of  the  cylinder.  Electron  energies,  monoenergetic  at  any  in- 
stant, change  continuously  throughout  the  pulse.  The  peak  density  of  the 

2 

emission  pulse  is  11.6  amp/cm  occurring  at  about  50  nsec,  while  the  peak 
emission  electron  energy  is  180  keV  occurring  at  about  6 nsec.  The  emission 
current  density  and  energy  are  such  that  substantial  space-charge- limiting 
would  be  expected  if  the  cylinder  were  evacuated;  instead,  air  at  200  pm 
pressure  provides  enough  ionization  under  these  conditions  to  effectively 
cancel  the  fields,  resulting  in  almost  no  deflection  of  electron  paths. 

Numerically,  the  problem  was  modeled  by  breaking  the  cylinder  into 
5 axial  and  5 radial  zones,  employing  a particle  time  step  of  0.2  nsec. 

The  code  calculates  a light  time  step  based  on  the  minimun  zone  size. 
Particles  were  emitted  from  15  radial  positions.  The  arbitrary  emission 
current  pulse-time  history  and  energy-versus-time  curves  were  modeled  with 
9 and  10  point  pairs,  respectively.  Mini-prints  were  requested  every  2 nsec 
and  spatial  prints  every  5 nsec.  Most  of  the  available  plot  options  for 
time  histories  were  requested. 

The  numerics  and  modeling  employed  here  should  not  be  taken  as 
representative,  especially  the  small  number  of  spatial  zones  and  the 
straight-ln  primary  electron  emission.  These  conditions  simply  provided 
a convenient  sample  problem. 

5.1  SAMPLE  PROBLEM  INPUT  CARDS 

The  input  cards  necessary  for  the  sample  problem  described  above  are 
shown  in  Figure  2.  Notice  the  comments  at  the  beginnings  of  the  data  cards. 
These  are  for  the  programmer's  convenience  and  are  not  required  by  the  code. 
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Figure  2.  Sample  problem  Inputs  for  cylinder  15  cm  long  by  30  cm 
diameter  at  pressure  of  2.63  x 10  '*  atm.  A time- 
dependent  energy  spectrum  Is  specified. 
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5.2  SAMPLE  PROBLEM  OUTPUT 


Selected  code  outputs  essential  for  interpreting  calculations  are 
given  in  this  section.  In  particular,  the  emission  particle  character- 
istics, the  mini-print,  and  the  spatial  print  are  described  fully. 

Definitions  of  the  plot  titles  are  also  given  in  the  glossary  in 
Section  5.3. 

The  sample  particle  emission  printout  is  given  in  Figure  3.  This 
printout  lists  information  pertaining  to  sample  particles  which  represent 
large  numbers  of  emission  electrons.  Initial  positions,  velocities,  etc., 
are  given.  In  cases  where  time-dependent  spectra  are  being  calculated, 
this  printout  can  be  specified  at  various  times  to  check  particle  emission 
against  input  quantities  [see  K0PT(3),  card  2].  Variable  names  are 
defined  in  the  glossary  (Table  4).  The  sample  mini-print  is  found  in  Figure  4. 

The  sample  spatial  print  is  shown  in  Figure  5 for  a time  of  5 nsec. 

Electric  fields  have  already  been  substantially  dissipated  at  this  time 
by  large  conductivity  currents  due  to  secondary  electron  migration  (drift). 

The  primary-electron  radial  current  density  is  zero  because  the  electrons 
were  Injected  straight  into  the  cavity  and  also  because  the  magnetic 
field  was  not  permitted  to  influence  electron  motion  [KOPT(37)  = 1]. 


23 


O x<  9 3 9 4 * 


--j'ooy'"''  — o*^>“f*'OOwnj^9 
I9*^0«SS'«OC'«^0'0-*<^ 
O ^ C f ^ 9 — * 

j*— ^OC-NJO** 

OO^OSOOO^OOOOOO 


^5999999999999999 

^^coococcc«co^o-o 


9 9 9 9 9 9 9 


'9  9 9 9 9 


999999999999999 

OOOOOOOOOOOOOOO 

♦ ♦♦♦♦♦♦♦♦♦♦♦*♦♦ 

a.»9-»»»»99‘99'9'a>9'» 


OOOOOOOOOOOOOOC 


9999999999999 

j>y>X'9'X'X'9'»/'j^X> 

CCCOOOCCOCCOf) 


9 9 
X /' 

o o 

9> 


9999999999999: 
^^0003  3000000C 


X"**  9 

9’'<%‘>jr«..'yf«-/\.  C-^0-*C^O  — 
Z9  r\i  »«^|\»'XJ»^|*09  9i/^knx 

0009X^CC9X‘\;-r09® 


ooooooo« 


' 9 9 9“  kT  .r 


O O ( 


XOOOOOOOOOOOCOOO 

^ocoo  oe  oooc.c.'^oeo 

oJ  9 OOOC-OC  OOOOOOGO 
XXXCOCXCXCOXXO-0 


> O O 9 O O c 


aOOOOOUCOOCiOUOOC 
^ OOOOOOOC  OOOOO  90 
^fvooooeoo  9000000 
i9'k/'J^J'kr9^9'J^9'‘AJ%J^kr'J'9' 


^ o o »•  - 


mOooooo: 

OCC  C OOOOC-OOOC  0<9 

• ♦♦♦♦♦♦♦♦♦♦•♦♦♦ 

90  OOOOOOOOOCO  90 

cooooooooooooooo 

99  9009090  OkT^^J^J' 
yk^(M-%9./X0f^  — 


•S40C  00900  9 0 000000 
OO  90090  0 0000  9C  O 
0099000090CCOOO 


9 ^ X9  O-VIXJ*^  : 


<u 

a 


X 

0) 


D 

o 


}-l 

CL 


CO 

CJ 


Li 

<D 

u 

o 

CO 

M 

CO 

4= 

o 

c 

o 

•H 

CO 

CO 

a; 

(U 

i-H 

u 


u 

cd 

CL 

t-J 

JH 

s 


0) 

M 

o 

00 


/ 


m 


Characteristics  are  given  for  a time  of  .2  nsec  or  at  end 
of  the  first  time  step. 
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Figure  5.  DYNACYL  spatial  print  e::ample.  This  printout  is  also  referred 
to  as  a "2-D  print"  elsewhere  in  this  report. 
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5.3  DYNACYL  VARIABLE  GLOSSARY 


A glossary  defining  DYNACYL  variables  and  printout  headings  is 
listed  in  this  section.  Between  this  glossary  and  the  Input  Cards  section 
given  previously,  all  of  the  code  input  variables  and  most  of  the  output 
variables  are  defined.  A few  output  variables,  primarily  of  a debug 
nature,  probably  remain  undefined,  but  all  essential  outputs  are  described. 
In  addition  to  the  "visible"  code  variables,  many  Internal  quantities  are 
also  found  in  the  glossary.  Thus  the  tabulation  is  useful  in  aiding  the 
programmer  to  understand  many  of  the  essential  calculatlonal  quantities  of 
the  code.  Many  variables  tabulated  will  be  found  to  pertain  only  to  the 
quasi-static,  or  Green's  function,  treatment  contained  in  DYNACYL  as  a 
subset.  The  programmer  should  keep  in  mind  these  are  simply  options  and 
in  no  way  indicate  DYNACYL  is  limited  to  quasi-static  situations. 
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Table  4 
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Sl.nSS*B»  OYNsm  ViMUHLtS  4^D  OUTPUT  MtAOINUS 
RADIUS  0^  CYLlNDtB  izH^  (Cw) 

angle  BET*EEN  YfLOCITY  Of  EPITTtO  ELECTRONS  AND  INCU“INr. 
PHOTONS  OR  the  surface  N0R«AL  DEPENDING  ON  KOPTIT)  (DEG,) 

AREA  OF  SIDE  Of  cylinder  (CH*»?) 

AREA  OF  TOP  OF  CYL  (CH**?) 

total  Particle  charge  in  cylinder  divided  cylinder 
vcilu’'e,  plasma  Charge  not  included,  (coul/crj) 

BACK  EMISSION  CURRENT  fraction,  E«PBESSED  AS  FRACTION 

Of  Current  emitted  frl'r  front  face 

RAGNETIC  field  at  cylinder  SIDE  aalL  AT  hDTTOR  (FORaaRd 
E“ISSION  FACE  J (GAUSS) 

tire  DERIVATE  Of  RAGNETIC  field,  (GUASS/NSEC ) , PROBABLY 
Ton  NURERICALLY  NOISY  TO  BE  OF  USE  IN  RANY  SCL  SITUATIONS, 

ERITTED  current  density  falls  off  by  factor  of  BFALL  by 
HBR^  IF  kOPT(B)*1 

PEAK  aRPLITUDE  RaGNETIC  field  ANYAHfRE  along  CYLINDER  SIDE 
aALL  at  this  Tl«f  (GAUSS) 

IBOTTCIR,  SIDE,  TOPI  DESIGNATES  The  NE»T  3 VARIABLES  IN  THE 
"INI-PRINT  list  ahIch  are  The  CURRENTS  striking  the  HOTTOR, 
SIDE,  AND  TOP  OF  The  CYLINDER,  RESPECTIVELY,  DIVIDED  BY 
THE  PEAK  ERISSION  CURRENT  FRO«  The  FORaARD  FACE,  USUALLY 

referred  to  AS  fractional  transmissions,  also  printed  Out 
labeled  fJBOT,  EJSIOF,  FJTOP,  The  Three  values  should  ado  to 
The  value  for  SURFJ  DEFINED  BELOa,  ( 0 I “E NS  I ONLE SS ) 

PrIRAry  electron  CmaoGE  IN  Each  CONE  (COuL), 

ALSO  USED  AS 

Charge  on  each  ERIttED  particle  in  TI“E  step  Print  or  in 
PARTICLE  description  print,  u/m  SARE  as  fob  electron  (CUUL) 

PRIRARY  electron  charge  density  (C0UL/“3) 

AHEN  successive  Tfrms  bECORE  less  Than  CNVRG  * total  Su“  of 
TERRS  UP  TO  That  Term,  Terminate  series 
FUR  GREEN'S  function,  ( 0 I RE NS I ONLE SS ) 

UUASI-STaTIC  ONL" 

PARTICLE  time  step  number 

APPRUA,  TiRf  Increment  for  pathes  of  particles  in  pictures, 

(SEE  DTF  ) EGUALS  DTF/NPLOT 
time  step  (NANOSECONDS) 

TIRE  STEP  IN  seconds  (SEC) 

DEL2( I )*E ( I»1 )./( I ) , III  TO  NZ  (C“) 

DR( I )lBII»I ).R( I ) , I>1  TO  NRRl  (CM) 

inject  charge  aT  DRANG  INTERVALS  IN  « ACROSS  TOP  Of  CVl,(CR) 
distance  of  source  from  cylinder  end  (CM) 

ELAPSED  t«ECUTION  TIRE  SINCE  LAST  CALL  To  GET  TOTAL 
EAECuTION  time  (MINUTES) 

fULLOa  parts  from  top  Uf  CYL,  FOR  PTF  NANOSECONDS 
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BEST  AVAIttBlE  CO?t 

Table  4 (cont.) 


light  TIhe  step  prom  T(N*i/^)  to  SEE  DT\2,(SEC) 

OTNi  light  TI"E  step  EHEIH  T(N)  TC  T(N*1),  T(1)iO,,  see  DTi».(3EC) 

OTSTEHd-VTSTEP)  SEE  TSTEP 

OxK»T  INTEBPOLiTE  TU  EINO  HELDS  FOR  ABS  ( ( z-z  0 ) /L ) , L T , 0 « I C » I T FOR 

hEThoDei,  abS((H-H0)/A)  ,LT,  DaICRT  FOR  “EThoOrJ 
If  IUPT(15)«0. 

EhaR  average  EfiERGV  OF  The  PRIhaRV  ELECTRHh  SPECTRU"  INJECTED  INTO 

the  CvL,  (hEv) 

EftARPT  ENER&v  of  pin  '"ITh  largest  MJhbe  R OF  privary  elEC'RONS 

being  evicted  at  This  tire  fro"  forward  face,  (“evi 

EBART  ERlSSinN  ELECIRUN  SPECIROH  AVERAGE  ENERGY  aT  FORhARO 

FACE  at  this  tire,  (wEv) 

EFRl  RADIAL  E*FIEL0  seen  by  FOLLCdED  PART,  AT  GIVEN  TJME  (V/R) 

EFZl  a«IAl  ErFIElO  seen  by  followed  part,  at  given  TI“t  (V/R) 

ELN  ELN( j )-ELN(NSPD)  are  The  NURBER/AREA  of  electrons  at 

ENERGIES  ENRGE  n*f  NRG(NSPO)  AT  ANGLES  PLUS  AND 
RINGS  ANG,  AT  RADIAL  PCSITIflN,  BANGED*  RANGENEPTS) 
ENljMBER/CR*«2)  ETOtAL  NuRRER/CR**i  FOR  A GIvEN  ENERGY  PULSE) 


EMISSION  POINT 
E“0 

energy 

ENRG 
E PSO 

EG 

ER 

E RMIN 

EVP 
f VB 

EE 

EEbAk 

EZRA* 

FCFSRV 

FCOLB 

fcolt 

F INJf » 

F JBOT 
F JSIDE 
F JTEIP 


I POSITION  OF  emitted  electrons  EC”) 

total  charge  E“ITTF0  into  The  CYL,  at  a given  radius  and  angle 
during  TmE  entire  pulse,  ecoul) 
initial  electron  energy  EREv), 
electron  energies  EREV) 

EPSILON  ZERO  EfrisSIVITy  of  VACUUM)  8,8SUE•1^  FAHAU/R 
FARAD  * COUL/VOLf 

SAME  AS  I EMITTED 

radial  electric  field  evolt/R) 

PEAK  amplitude  Radial  electric  field  anynneRE  along  cylinder 
side  *all  at  This  tire  evoli/r) 

z velocity  of  emitted  particles  er/seo 

Radial  speed  of  fritted  particles  Em, SEC) 

AFIAL  ELECTRIC  FIFLD  EVOLT/m) 

‘»I*L  electric  field  On  aajS  AT  foRpARd  FRISSON  fACEEVOLT/m) 

a»Ial  electric  field  on  a»IS  at  top  face  Eapay  fRO“  foRpard 
EMISSION  FACE)  EvOlT/m) 

fractional  NON.CE'nSEBVaTION  CF  energy  EO  fob  complete 
CONSERVATION)  EDImENSIONLESS) 

LIRE  FCOlT  only  for  bottom 

calculate  FRAC  ThaNS  Thru  top  for  R ,LT,  and  ,GT,  FCOLTfR^ 
emit  electrons  beTpeen  »»0  and  RfFINJEp  • R^ 

SEE  BST 
SEF  hST 
SEE  BST 
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BESr  AVAILABLE  COPY 

Table  4 (cont.) 


^OHOT 

MJSiDfc 

► otop 

f T 

ftbut 

FTSiDfc 
FT  TUP 
F»*FLD 

Ft*.FLD8 

F«FLPt 

FwIN 

UNGPD 


H«Pn  I 


FRACM'IN  UF  total  primary  ELFCTHflN  CHARGE  *nhICh  HAS  STRUCK 
CVLINOfH  BDlTO"  F4Ct  UP  TO  I^IS  TI"t  ( D I “t NS  1 O^L^ SS ) 

Fraction  Of  total  pbimjhy  electwon  cbarce  «hich  mas  stboc" 
CyLINOEP  SlOf  y.ALL  up  to  this  time  (DIMENSIONLESS) 

FBACTION  Of  TOTAL  PBIMABV  ELECTRON  CHARGE  hMKM  mAS  STBliC" 
cylinder  top  face  up  to  THIS  time  (DIminSIunlESS) , 

relative  emission  current  pulse  hEIGTH  at  present  time 
(NORMALIZED  TO  peak  value  OF  (,)  FROM  THE  EORkARO  EMISSION 
FACE. 

fraction  of  energy  injected  up  to  This  Tjme  mhICh  has  left 
Thru  top,  (DIMENSIONLESS) 

like  FTbUT  but  for  side 

LIKE  FTBOT  but  Ff!R  TOP 

fraction  of  total  energy  inside  CYL,  *mIch  is  in  fields 
(DIMENSIONLESS) 

fraction  of  field  energy  in  b-field  (DI“ensionless) 

fraction  of  field  ENeB&y  in  E'FIELD  (DIMENSIONLESS) 

fraction  of  total  energy  injected  up  VQ  this  time  which  is 
INSIDE  The  CYL,  AT  THIS  TIME,  (DIMENSIONLESS) 

Radial  and  a«ial  fields  due  to  annuli  or  rings  of  unit 
Charge  depending  on  I0PT(15)  (STAT.VOLT/Cm/STAT-COUL  f 
R,E»t3  VOLT/M/COUL) 

FLAG  FOB  ANG,  DIST,  SET  FROM  NANG  VALUE 
1 ANG,  DIST,  ON  CARDS 

0 cos,  ANG.  DIST 

-I  ISOTROP,  ANG,  DIST 

MAGNETIC  field  (AMP/m) 


I emission  index 

I EMITTED  total  primary  electron  CURRENT  EMITTED  AT  THIS  TIME 

Includes  both  emission  faces,  (ampsi 


IBOT 

IBOTIN 

ICHECK 

IDL 

IDT 

IDYN 
IF  IN 

INJEk 


total  primary  electron  current  striking  cylinder  bottom  at 
This  time  (AMP) 

total  primary  electron  current  striking  bottom  face  oe 
cylinder  at  radius  less  than  FCOlB  * cylinder  radius. 

FCOLB  SET  AT  beginning  OE  SuEIROllTINE  TOIMRZ,  CURRENTLY 
EOUALS  ,(i7,  (AMPS) 

CHECK  convergence  of  SERIES  EVERY  IChECK  Terms  (SEE  IPHINT) 

IF  ICMECF  input  as  negative,  check  convergence  of  sums 
BY  SUMMING  IN  reversed  order,  (IOPT(J)) 

increment  in  The  INDEX  OE  The  PARTICLES  WHICH  WILL  BE  FOLLOWED, 
There  will  be  NSB0*NEPTS«2«NANG  particles  emitted  each  time 
STEP  wITh  angle  varying  EIBST.TmEN  RADIOS  HE  EMISSION,  Then 
energy,  see  koPT(Z) 

INDEX  increment  In  Time  for  storing  fractional  transmissions 
AND  potential  VS,  time  (SET  TO  NT/NJDT  « 1) 

0 DoASIwSTATIC  greens  FCN  method 

1 Dynamic  max.Eq,  method 

1 IF  ,lT,  10  SEC  XQT  TIME  remain  OR  IE  ON  L*ST  TIME  STEP 

0 otherwise 

total  number  oe  particles  injected  into  cylinder  op  to 
this  time 
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BtSrAVAW^ 

Table  4 (cent.) 


IOPT(0?) 


StT  fOU«t  TO  IPRINT 


lUPT(OJ) 
lOPT (OU) 

lOPT(OS) 

lOPT (06) 

lOPT (QT) 

lOPTfOS) 
!OPT (OR) 

lUPT ( 1 0) 
lOPT (11) 
lOPT ( 1?) 

inPT(lJ) 

lOPT ( lu) 
ItiPT(lS) 

I0PT(16) 


SET  E(JU^t  TO  ICREC« 

Plot  sums  in  series  e»R4nsiun  vs,  n eor  single  sur  sols, 

FOR  EVERY  lOPT(u)  IR  UNO  12  viLUtS 

PRIM  OUT  every  eCN  eor  E^CR  n for  single  SU"  SERIES 
FUR  Every  iopt(5)  ir  and  iz  values  (see  ioptid) 

,lt,  0 print  only  sum  vs  n eor  Every  iars  (iopt(s))  ir.iz 

.GT,  0 REA[1  in  I0PT(6)  (ZO.RO)  pairs  (“UST  be  values  included 

AS  FIELD  PTS.)  AND  CALL  FIELDS  FUR  THESE  (ZO, RO)  ONLY 

PLOT  EZ.ER  vs  Z and  R for  every  I0PT(7)  ZO.RO  PAIRS 
Putting  ioptib)  curves  on  each  graph  (if  iopt(7),gt,o) 

Put  niPT(8)  CURVES  ON  each  graph  IE  I0PT(7),GT,  0 (max  of  u ' 

,GT,0  RAP  EW.eZ  INTO  33  INTERVALS  (3-D  PLOT) 

,LT,0  rap  ER.EZ,  ''Lx  no  of  ITERATIONS  INTO  J3 
intervals  (3-d  PLOT) 

,gt,  o take  Divergence  of  e fiEld  everywhere  and  print  out 
fur  each  20, ho  using  Subroutine  divErg  (for  individual  rings) 

,GT,  0 PRINT  OUT  eP.EZ  as  FCnS  OF  «,Z  FOR  EVeRT  lOPT(ll) 

RO.ZO  POINTS 

,&T,  0 integrate  greens  functions  fur  fields  over  charge 
distribution  which  Is  specified  in  etiital  for  debug 
Purposes  only 

,NE,  0 call  ChErE  to  calc  sums  for  ER(TOTaC)  and  E2(T0TAl) 
wJTh  constant  RhO  assured 

.LT,  0 take  divergence  of  fields  CALCED  with  chEkE 
.LT,  -1  EXIT  after  call  CHEKE 

.GT,  0 PRINT  OUT  greens  FCNS  FOR  ER,EZ  FOR  EACH  R,Z,RU,ZO 

0 CALC  GREENS  FCNS  FOR  INFINTESIRAL  RINGS  OF  CHARGE 

1 CALC  FIELDS  DuE  TO  FINITE  vOLuRE  ELEMENT  OF  CHARGE  FOR 
GREENS  FCNS,  IE  SMEAR  OuT  THE  CHARGE  OvER  A vOLURE 
ELEMENT,  GOOD  FOR  methODkS  Only 

,CT,  0 FAST  DRUM  UNIT  FOR  GREENS  FCN 


U)PT(17) 

.GT,  0 

Kt  AO 

WtSTAHT 

TAPt  ON  UMT 

lOPTdT) 

lOPTdS) 

.GT,  0 

Wf.STAPT 

TAPe  ON  UMT 

lOBTdS) 

IPLOT  Z print  particle  positions  and  velocities  vs  time  for 

SELECTED  particles  , AlSO  PLOT 
1 Plot  particle  paThes 

IPRINT  IF  ,GT,  0,  intermediate  print  of  SERIES  EVE»Y  IChECk  TERRS 

IF  ,GT,  -1,  PRINT  NO,  OF  TERRS  ySED  FOR  EACH  R,  2 
CALCFD,  (IUPT(2)) 

IPSTOP  NUMBER  OF  PxRTICLES  TRjNSFER»ED  TO  CONTINUUM  PLASMA  DURING 

This  time  step 


iHMjx  Parameter,  no,  of  rxoial  zones  ♦ i 

iRUMAx  Parameter,  no,  op  rxoial  zones. 


ISIDE  total  primary  electron  current  striking  CTlINDER  side  at 

This  time  (AMP) 

istopt  cumulative  number  of  Particles  transferred  to  plaS"a  up 

TU  this  time, 

ISTURE  tdIM  0 DONT  store  part  positions, etc  this  time  step 

1 STORE  PART  POSITIONS, ETC  THIS  TIME  STEP 


i ' 
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ITJme 
I TOP 

ITOPIN 

IZWftX 

J 

J tMITTtO 
JtiUT 

JOT 

JNOa 

JOUTR 

JOUTS 

JOUTT 

JPt  AH 

JR 

JR  conductivity 

JSIOf 

JTOP 

J^ 

j;  conductivity 

WPFLD 

KPULSfe 

KH 

KZ 

I 

u 

LL 

LOvvep  aall 
CHARGfc  density 

tPt 

LPm] 


COPY 


Table  4 (cent.) 


PSBTICLt  TIMf  ST(P  MlWHfK 

TuT4t.  PB1“A«»  ^LtCT»n^  CUPPE'iT  STPIKINO  C’LINDEB  tup  at 
This  time  (AMP) 

Same  as  IBQTiN  Of-LT  FOH  CTLINOEP  TCP  (FACE  away  F«0M  EUPaABO 
EMISSIUN  FACE ) (AMP, 

PABameTEB,  no,  (IF  axial  Z(m£S  ♦ 1 

INCE* 

TUTaL  CUBBENT  at  TIMf,T,  EMEBING  CYL,  (AMP)  SEE  EO 

IboT  DIYIOEO  HY  bottom  FACE  AREA,  GIVES  Average  pbimapv 
ELECtBON  Cu®Bt^T  density  incident  On  bottom  face,  (AMP/CM^) 

Running  incbe“Ent  in  time  fob  storing  fwactional  transmissions 
AND  potential  vS.  time 

EMISSION  electron  current  density  at  axis  of  forward 
EMISSION  face  at  This  time  (AMP/CM2J 

Same  as  JDuTt  hut  for  bottom  face 


same  as  JOUTT  But  for  side  wall, 

PRIMARY  electron  current  DENSITY  STRIKING  CYLINOE*  TOP 

Surface  in  eacb  Radial  zone,  (amp/cm^. 

Peak  emission  ELECTRON  CubbE'iT  density  at  axis  Of  forward 
EMISSION  face  (AMP/dm,;) 

primary  electrcin  Radial  current  density  (A“P/Ma) 

Radial  current  density  in  each  zone  due  to  Secondary 
electron  migration,  (AmP/m^j 

ISIDE  divided  by  cylinder  SIDE  "ALL  A«EA,  GIVES  AVERAGE 
primary  electron  current  density  incident  on  Cylinder 
SIDE  WALL  (AMP/C'V) 

ITOP  DIVIDED  BY  cylinder  Top  face  area,  GIVES  avERAgE 
primary  electron  Current  density  incident  on  top  face  of 
cylinder  (ampycmZ). 

PRIMARY  electron  axial  current  Density  (amp/mJ) 

AxIAL  CURRENT  density  in  Each  zone  due  to  secondary 
electron  migration  (AMP/m^i 

PRINT  flag  used  in  TdIM,  0 NOT  A PRINT  CYCLE,  I P«INT  CYCLE 

FLAG  used  to  indicate  max,  storage  for  no,  of  parts,  has  been 
reached.  p Turn  pulSE  OFF,  i lEavf  pulSE  QN, 

Surface  current  density  in  radial  direction  (amp/m, 

surface  current  density  in  axial  direction  (a“p/m) 

length  of  ctlINoER  L«ZZ  (CM) 

parameter,  max,  no,  oE  particles  which  cak  “E  I.JCC.ti  AT 
EACH  injection  time,  (?«NANG«NtPTS.NSPU,LE ,LE ) 

Parameter,  max,  no,  of  particles  code  can  follow. 

Surface  charge  density  in  each  zone  on  lo„er  Cylinder 

FACE  (C('UL/MZ) 

parameter,  no,  of  axial  Zones  ♦ i 

PABTICI.E  initial  angular  momentum  in  The  unfortunate 
units  of  G-CmwmysIC  whICh  means  The  number  is  I0«*s  higher 
T"E  mkS  unit  of  JOUlE'SEC, 
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IU< 

LKG 

L^»AT^l-^PG) 

makks 

ThOO 


MJl;- 

MM£,h.S 

^A^G 


dicta  W is  W^/Sijw  Cf  C»ATIU(l)  Thru  C'^A  T 1 0 { ^hCP  ) , TnfcN 
DEL«{1)  « L«aTIO(I)  • »*I\,  DfecTA  « 

pARA'^eTfeR,  f.n,  HP  Radial  zc\es  ♦ \ 

5A«t  AS  LRaTjo  t»CfRT  pqr  axial  zones 

^A*,  NP,  Of  TfeR^'S  ALUO^ID  tAC^^  SU“  IN  Tnt  StWItS  POP 
GRtfe^'S  PuNCTION,  SFT  TU  -RAiNS 

I IP  L/A,&t,i,  SINGLE  SU''  SPL,  -ITh  SI'nh,  CUSm  FCNS, 

Z DUuP'LE  Sjp*  SOL,  ^IT»-  OttniNAHf  h^SSEL  pCNS,  DE?1U& 

CPTIHk,  vuST  be  Sf T TN  PCRTRak, 

! IP  L/A,iT,i,  SINGLE  Su'^  SCL,  *ITh  MpoiPlED  bESStL  PCNS, 

• I DO  single  Sl>R  8PLS,  PO«  “ETmOO»1  and  5,  OE«U& 

CMTI^3^,  must  he  set  in  FORTRAN, 

PARA’^tTE®  variable,  EmUALS  NJUT 

paRa-'EtE®,  pp,  p TE‘»''S  alL  VEO  in  series  H)R  greens  PCn, 

NO.  pp  angles  a»  Particles  are  e'^itted,  see  le 

read  i»  In  ■pPTdO) 

,GT,  0 CODE  reads  AN&,  »)IST  f’N  CAwr:  aITh  ^A^G  OiNS 
,lT,  0 coot  sets  cos  ang,  OIST  aITh  iahS(Nang)  bins, 

,lT,  .100  CODE  set  ISOTROP,  ang  OIST  r,lTH  lAfaS<NAN&)  .100 

angle  hins 

»,z  r.EOw,  • NO,  pp  Radii  at  #.mich  incoming  electrons  are 

ERJTTfc  ■ FHUR  The  top  OP  TnE  C^L,  (SEE  LE) 
read  in  in  «PPT{i<i) 

input  ’m(  incident  electron  SPEcTRuR  evert  NE3«1P  time  steps 
(RE  ad  in  in  kOPT  <a) ) 


net  C^'ARGE  net  C'’AR:E  uENSdr  IN  tAC«  ZONE,  including  primary  and 

density  secondary  electrons  and  ions,  (COUL/mJ) 


NPOLPyi 


Parameter  varUhle,  no,  qp  particles  TO  he  poll^'^Eu  pgr 
trajectory  Plot 


MNJEN 

NJDT 

NLF  AVE 

NL  vTPT 

NNR 

NNWQ 

NNZ 

NNZU 

NO , LE  A V 1 NG 


NO, UP  angles 
NO,  tip  particles 

NO,  OP  PARTICLES 
IN  each  zone 

NPaRT 


total  number  pp  Particles  injected  into  ctlinoer  up  to 
This  time,  presently  correct  ONLY  IN  forward  emission  cases, 


STORE  and  Plot  practiunal  tRanshissIons  and  potentials  po»  njot 
time  values  MET^^EEN  0 AND  TRA*,  see  ICT,  SET  TO  mjoT 


tuTal  np,  op  Parts,  leaving  c^'l,  P-  time  increment,  delT,  at 

time,  T, 


total  number  op  Particles  striding  ctlinder  •'alls  up  to 

This  time. 


PARADE  ter 


^0,  pp  PIELO  PT5,  IN  radial  direction. 


PARAME  ter 
PARAME  TE  w 

parameter 


NNR 

no,  op  field  ptS,  in  axial  direction, 

NNZ 


TotAL  Number  op  Particles  leaving  system  by  striking 
cylinder  •'ALLS  during  PREVIOUS  Time  step,  does  nut 
INCLUDE  particles  converted  TQ  PLASma  BACKGROUND  IN  AIR 
pressure  cases. 


NO, OP  angles  at  r.HlCN  ELECTRONS  ARE  INJECTED 

Number  pp  particles  being  POlLOk-ED  at  present  time 

number  up  Particles  representing  p«i“apy  electrons  in  each 

SPATIAL  ZONE,  DOCS  NOT  COyNT  PARTICLES  OuTSlDE  CYLINDER 
but  aIThIN  l/^  ZONE  hhICh  aRE  BEING  PULLn^ED  FOR  PROPER 
boundary  treatments, 

nu“be«  op  Particles  being  polloreo  this  time  step, 

ALSO  printed  Out  a8  LAST  NuRBER  on  top  line  of  TmE  MIM.PRINT 
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Table  4 (cont.) 


^PG 

NO,  Of  iXUL  ZQNtS 

NPGPl 

♦ 1 

NPHJ 

Stt  KOPT(IO) 

NPl.0^ 

PaBAMjTfR  V4P]4WLE  4hICm  MUST  ,C.l,  NO.  OF  POINTS  MtlNC, 

FLClTTtD,  I.t,  NPUOT  ,St,  M4»(NP,NZ,M4XNS)  IN  tflELD, 

In  TDIm,  NPLOT  Is  4 p4M4MHt«  VARUetf  OlidNl.  THE  NO,  OF 
PARTlCLt  POSITIONS  TO  BE  PLOTTED  for  OTF  NANOSECONDS, 

NPLOT2 

P4R4METEB  V4RI4HL!  E0U4L  TO  NPLOT«2 

NR 

NO,  OF  FIELD  PTS  IN  R4DI4L  CI»ECTION 

NRAOI 

FL4G  USED  TO  signal  more  INCOMING  electrons  4T  anOThEW  angle 

0 NO  MOPE  electrons 

,NF,  0 MORE  electrons  at  ANOTmEM  ANGLE 

10  MriNOCMROMaT  IC  BEAM  aITh  naNGaNEPTS*?  PARTICLES  INRuT 

Into  tor  of  Cyl,  EACm  step 

NRG 

NRGR* 1 

NfiGt* 

NO.  OF  RADIAL  zones 

NRGPPl 

NKGR‘1  EDUAlS  NRG 

NRMJ 

NR-l 

NRO 

NO,  OF  R POSITIONS  OF  RINGS  OF  CHARGE  BEING  CALLED  ( ,GE,  1 

NRl 

number  of  particles  reflected  Through  c’lInder  axis  this 
time  step 

NR? 

number  of  particles  hhICh  STRUC*  CYLINDER  SIDE  HALL  THIS 

TIME  STEP 

NSPO 

NO,  OF  ENERGIES  IN  spectrum  of  INCOMING  ELECTRONS  (SEE  LE ) 

NT 

NO,  OF  time  steps  to  be  ta«En,  nT»Tmai/0ELT 

NTOT 

Total  no,  of  electrons  entering  thru  top  of  cTl,  In  a GIyEN 
IF  ZERO,  calculated  from  jpear, pulse  area,  and  Cylinder  area 
EuuALS  total  no,  injected  in  hhole  Pulse  internally 

BUT  INPUT  AS  NQ,/CM.»? 

ntpuls 

NO,  OF  time  and  pel,  pulse  height  VALUES  USED  TCI  DEFINE 
EMITTED  Current  pulse  (roptcud 

NTSk IP 

PRINT  OUT  EVERY  ntSktr  Tjmf  STEPS 

READ  IN  IN  KOPT(IS) 

N T S T E.  P 

NO,  OF  DIFFERENT  time  STEP  SIZES  (READ  IN  IN  aORTCIZ)) 

N? 

NO,  Of  FIELD  PTS  IN  axial  dIRECTiov 

NZMi 

NZ  • 1 

N/0 

NO,  OF  z POSITICFS  OF  RINGS  OF  CHARGE  FOR  GREEN'S  FUNCTION 

NZ| 

number  of  Particles  hhich  struc  cylinder  bottom  this 
time  step 

nZ? 

NUMBER  OE  PARTIClES  nhICh  sTRuC"  CYlInDER  top  this  time 

STEP 

OUTfcR  A all 

Surface  charge  density  in  each  zone  along  CTlINdeR  side 

CHARGfc  density 

hall  (CUuL/Mj>) 

PA 

AIR  PRESSURE  IN  cylinder,  ALSO  FLAGS  UUASI*ST4TIc  CALCULATIO 
IF  SET  TO  •1,E*Z0  and  RROHlE"  must  be  vacuum  only. 

PHI 

particle  velocity  initial  azImuthal  angle,  phIio  means 
particle  initially  goes  haDIally  OoThARD,  PhI«ipo 

MEANS  particle  initially  COES  INaARD, 
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Table  4 (cont.) 


PLASMA  OSCIlLA. 
-7IUN 


POTL^MAL 
PEL,  *r'sA 

«eL,  ^0  ZSL 

PWl'^Afrv  fcLECTPO^ 
lONWATjON  PATe 


APPPOXlMftTE  tXPPkSSIO^  FOP  PlaS^-A  OSCILLaTJONS  HASEL) 
AVfPAGE  ChAPGE  OE^SITr  IN  C^Ll^OEPt  SmOmLO  HE  VA&UfctV 
consistent  isITh  (IhSEPVLD  LO^  FPEHUENC''  PI^G1NG  PfPlUO  IN 
Gu4SI*STAtIC  SITuATIHNS  (NSEC) 

ELECTPIC  potential  HPTilNED  pv  INTtGPATiNG  »A0IAL  ELECTRIC 
FlfLO  IN  FPOM  C'^Ll-vOEP  SIDE  •'■ALL  (VOLTS) 

ELfCTPIC  PfJTENTUL  OHTaINEo  RV  INrtG»ATING  AXIAL  ELECTHIC 
FIELD  DOisN  FPDm  TOP  FACE  OF  CruINDEP  (E»L)  (vOlTS) 

PATE  AT  PPI^APV  electrons  IONJZE  ThE  Altf  MOLECULES, 

(tLECTP0NS/«5/SEC) 


PULSE  (1*NTPULS)  PELATlVE  VALUES  OF  EMISSION  CUPPENT  DENSiTv  AT  TI«ES  TPuLSEd 

ThHU  ntpulS) 


0 

QDIF 


UE 

GJNJEN 


OINStD 


QLE  a vE 


GlEavp 


cmapge  on  Particles  e'^itted  this  tip^  step  (coud 

fractional  difference  het«^een  cu'^ulative  emitted,  and 
absuRBEO  PrIrAPY  electron  CHARGES  anO  The  Amount  inside  TmE 
Cylinder  according  td  the  ccoe,  value  should  be  s^^all 
(LT,1)  indicating  good  conservation  op  something  has  gone 
aPONG  aITh  the  calculation,  see  note  on  '•JINSyD 
(OIhfnSIONlESS)  , 

electron  charge  (CDUL) 

total  charge  t^’tTTED  l»P  to  This  TI‘*f  FRDM  ^^Th  E-'ISSIUN 
Surfaces,  (COul), 

total  primary  electron  Charge  •'IThin  the  c’fv  inoer  at  this 
Tl^'t  (COUL),  DOES  NOT  count  TmOSE  PARTICLES  f-nOSE  CENTERS 
ARE  Outside  the  cylinder  but  .sitnin  \/i  jone  hhic"  are 
fullO'^ed  for  obtainug  phopEr  rtOl'NpARv  valuEs  OF  Currents, 

total  Charge  -^hIch  MftS  struck  all  three  C'^lINdeR  surfaces 

UP  to  This  TIRE  (COUL)  INCLUDES  TnfiSE  PARTICLES  OUTSIDE 
The  CTLINOER  y-hICh  are  still  being  FOLLUaEO  to  \fl  ZONE 

AAAY  FOR  PROPER  tfOUNOARY  TREATMENT, 

Nt  T CALCULATED  PROPEPL^'  AT  PRESENT, 


qouTp 
DOUTS 
DOUT  T 


SAME  AS  OOuTt  hut  for  hOTTom  face, 

Sa*^E  as  QOUTT  muT  for  SIDE  •^ALL, 

total  primary  elec^^on  Charge  per  unit  area  *hic-  -as 
sTRyc*  Cylinder  top  surface  in  each  radial  zone  (Coul/chz) 


GR  1 


CUMULATIVE  amount  of  ChAR&E  REFLECTED  ThRQuGh  C'^LINDER  AXIS 

LP  TO  This  TlMf  STEP  (COUL) 


QR^ 

ORZP 

USm 

USTCiP 

DSTC'PT 

0/1 

(SZlF- 

0/? 


Cumulative  Charge  y-hjch  has  STRuc*  cylinder  side  hall  1.  R 
Tu  This  time  (COul) 

Same  as  nZ2P  tiUT  for  C'^LINdER  side  i^ALL. 

ChArge/mASS  for  electron  (CCUL/Gm) 

amount  uf  charge  converted  T(i  The  Plasma  formulation 
elfctuon  charge  density  During  this  time  step  (COud, 

cumulative  amount  of  CHAR(iE  converted  to  The  Plasma 
Formulation  electron  chahge  density  up  Tp  This  ti-e  (col'D, 

Cumulative  Charge  ahich  has  STruC«  CvlInoer  hOTTcm  yp  to 
This  time  (COUL) 

Same  AS  QZ^P  H-  T FOP  cylinder  bottom  face. 

Cumulative  charge  ahich  has  struck  Cylinder  top  u**  to  this 
time  (COUL) 


OZZP  Sum  of  ODuTT  (I)  TIveS  AREA  OF  RADIAl  ZONE  1 UvER  anOlE  FACE, 

PROVIDES  internal  consistency  ChEC"  OF  Cumulative  Charge 
Striking  top  Su‘*f*CE  ahen  compared  hIth  qz^,  (cuuD 


ito 


R 

H SPEED 

RADIUS 

RA^'C^^EPTS) 

RBAR  ( 1 ) 

RDD 
RE  FC 
RCi 

RGPt l»NRGP) 

RHO 

RhO.E 

RmO"! 

RINJEK 

RU( l^NWO) 

RPOS 

PO 

R1 

SECu^fARY  ELEC» 
TRON  U1NIZATION 
RATE 

SIOEMT 

SIGPA.M 

SIGMA.Z 

SPO 

SPEED 

STEP  NO, 

SUMF  J 

SUMF  0 


BEST  AVAILABLE  COPY 

Table  4 (cent.) 


PABTICLt  initial  madial  PCISITIUN  (CM)  IN  EMISSU'N 
CHAHACTtBISTICS  PHINTOUT  , ALSO  USED  fOB  HAOIAL  (.BIU  POINTS 
IN  SBEEMS  ECN  calculation, 

BjoIal  speed  up  emitted  ELECTWUNS  (M/SEC) 

BadIAL  position  OE  emitted  ELECTBDNS  (CM) 

BiNGd)*  BaNG(NEPTS)  aBE  The  BaDII  e BQM  ahIch  eLECTmoNS  ABE 
emitted  EBOm  The  top  OE  The  CVL,  (tVENL*  SPACED)  (CM), 

SEE  BINJEK, 

BEIAH(I)  S .Sa(B(I)a8(I41))  , Irl  TO  NBMl  (CM) 

ACCEL  IN  B direction  ( M/SEC**i) 

BEFLECT  befo  OE  the  ChaBGE  striking  the  CtL,  walls 
Radial  eone  boundaries  (cm) 

ZONE  center  BADII  (CM) 
charge  density  (COUL/CMAA3) 

BHOdZ.IB)  is  fob  zone  bounded  by  ZdZ)-Z(IZ*l),MdB)-BdB*l) 

SECONDARY  electron  Nijmeie  r DENSITY  (ELEC  TBOWS/M3  ) 

ION  number  density  dONS/M3) 

OUTER  RADIUS  OE  AREA  OVE»  WHICH  CHARGE  IS  INJECTED  INTO  CYL, 
(CM) 

R CUOROS,  OF  RINGS  OF  charge  IN  SHEEMS  FCN,  (CM) 

R POSITION  Of  FOLLOWED  PART.  AT  (,IVEN  TIME  (CM) 

FUR  lOPTdS)*),  RO  IS  inner  R value  of  The  finite  volume 
element  Of  Charge  of  ThE  source  annulus  in  The  GREEN'S  FCN, 

INNER  RADIUS  OF  SMALLEST  RADIAL  ZONE  (CM) 

RADIUS  OF  cylinder  (CM) 

RATE  AT  WHICH  sEConDART  electrons  IqnIZE  the  air  molecules, 
(ELECTRONS/M5/SEC) 


LIKE  RAKEMT  Only  for  SIDE 

RADIAL  CONDUCTIVITY  OF  PLASMA  IN  EACH  7UNF  (MhO/m) 

A»IAL  CONDUCTIVITY  OF  plasma  I w EACH  ZONE  (MhO/M) 

SPEED  OF  emitted  PARTICLES  (M/SEC) 

SPEED  DF  emitted  electrons  (m/SEC) 

PARTICLE  time  step  number 

SUM  OF  The  currents  SThIking  all  Three  cylinder  surfaces  at 
This  time  divided  by  the  peak  emission  current  from  the 
forward  face,  (DIMENSIONLESS) 

appru»imate  charge  conservation  indicator.  Should  be  close 
to  unity,  equals  sum  Of  charges  striking  Cylinder  bottom, 
SIDE.  AND  TOP  PLUS  Charge  inside  all  divided  by  total 
CHARGE  emitted  up  to  This  time,  does  not  count  charge 
OUTSIDE  cylinder  But  wIT«In  l/Z  ZONE  SO  IS  USuAllT  SlIQhTly 
LESS  Than  1,  Should  add  to  unity  when  emission  Pulse  is  over 
AND  ALL  primary  electron  CHARGE  MAS  LEFT  The  CYLINDER, 
unless  particles  Of  CHARGE  wERE  CONVERTED  TO  the  BACKGROUND 
PLASMA,  (DIMENSIONLESS) 
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Table  4 (cent.) 


T 


tl*'t  (\StC) 


Tl 

Tf  e 

Tuitn* 

TKfTi 

TIM[ 

IJMt  C.»ID 
T INSlfl) 
TmAK 


fhjctiun  Of  THf  total  ffissic'N  cu^wi^t  pulse  apea  uccopied 

Hy  Tt-E  PBESEMT  f.vlSSIO^  PULSf  xfT&TM  AT  TmE  FUHyAHD  FACE 

T1*'ES  the  particle  Tl^E  STER,  (DI“E^SIO^LESS) 

5A"E  AS  Te.  ONLT  for  back  E“ISSIUN  face, 

'"ExT  time  at  yiHjcM  Tn  pic«  u®  Parts,  fur  RLnTTiKG  patbes  (ns) 
Particle  velt'citt  initial  rclar  angle  reasurei)  frub  surface 

►URRAL.  (DEGREES) 

TIF'E  (NSEC) 

time  step  (NANO. SEC) 

total  kinetic  ENERGT  DF  eLECTRI'NS  inside  C’L,  (JOUL) 

VAX  TjuE  (NANli.SEC) 


total  charges 

On  The  STSTfv 


INSIDE 

Cylinder 


UPPER  BALL 


The  five  variables  LISTED  HELO*  PERTAIN  To  PRIMARY  ELECTk(1N 
PARTICLES  AND  FIELDS  ONLY,  The  NUMbERS  ARE  OBTAINED  From  Tmx 

CuNTiNuiTy  Equation  and  they  allo*  a check  dn  the 
conservation  of  Charge  in  the  portion  of  the  code  bhiC" 
CONVERTS  charged  particles  to  currents  ON  THE  SPATIAL  GRID, 
plasma  charges  and  charge  motion  are  not  included# 

TOTAL  primary  electron  charce  inside  Cylinder  (Coud, 

MINOR  inconsistency  RETbEEN  this  VALUE  AND  OINSYD  VARIABLE 
probably  due  to  finite  time  stop  SIEE  used  in  integration 
OF  continuity  eouation.  This  number  is  obtained  from 
Summing  array  for  charge  density  in  T"E  “AxbELL'S  EUUATIUN 
portion  of  TmE  code.  bmEreas  UINSYD  is  obtained  directly 
from  particle  storage  Arrays, 

NET  Charge  on  cylinder  top  face  (coud 


LObER  ball 

Outer  ball 
SYSTEM  total 


TPULSE 

TSTABI 


TSTep( 1-NTSTeP) 

UPPER  ball 
Charge  i.ensity 

VMaG 

VMAX 


net  charge  on  cylinder  bottom  Face  (coud 

net  charge  on  cylinder  side  ball  (COUL) 

SUM  OF  The  charges  inside  and  on  The  cylinder  balls.  Should 
ADD  TO  small  value  indicating  charge  conservation,  (COUL) 

times  at  bhIch  emission  current  pulse  is  defined  (NSEC) 

IF  reading  a dump  Tape,  RESTART  at  time  ,ge,  TSTaRT 
(CAN  BE  0,)  SEE  lOPT(lT)  (NS) 

time  step  is  DTSTep(I)  BeTbeeN  TsTeP(I)  and  TSTep(I«1)  (NSEC) 

SUPFACE  CHARGE  density  I Fi  each  EONE  (in  UPPER  CYLINDER 
FACE  (COUL/M?) 

MAGNITUDE  OF  Initial  particle  velocity  (m/seC) 

largest  magnitude  of  The  electric  potential  ANYbHERE  ALI'NG 
THE  Cylinder  axis  *t  this  time,  ( vOlts) 


vpmI 

VR 

VRI 

VEl 

VZ 

BF 


/>ARTICLE  initial  azimuthal  VELOCITY  (m/sec) 

PARTICLE  initial  b»0I*L  velocity  (m/sec) 

RADIAL  SPEED  OF  FOLLOWED  PaRt,  At  GIvEN  ylME  (m/seC) 
axial  speed  of  foLLObED  part,  at  GIVEN  TIME  (m/SEC) 
PARiicLE  initial  a«Ial  velocity  (M/sec) 
energy  in  e and  m fields  (jnuL) 


1*2 


»kI»yC(iw 


v^lNSVD 

Z 

2 SPftD 

ZHPtT 

zon 

2tI(l-^ZU) 

ZPOS 

20 

2\ 

2? 
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Table  4 (cont.) 

coBHtCT  tot»l  ENt»c.Y  «HicH  Should  he  in  cyl,  (Euoals  tNencY 
INJECTED  • energy  hmICm  mas  LEFT  UP  TO  This  TImE)  (JUUL) 

total  kinetic  ♦ field  ENER&y  inside  CVL.  (JOUL) 

AxIAL  ZONE  BOUNDARIES  IC")  . ALSO, 

PARTICLE  INITIAL  AXIAL  POSITION  (CM)  IN  EMISSION 

characteristics  print, also  grid  points  in  green's  fcn 

AXIAL  SPEED  OF  emitted  ELECTRONS  (“/SEC) 

ZBARd)  1 .S*(Z(I)TZ(I*I))  , Ixl  TO  NZMl  (CM) 

Plot  magnetic  field  at  outside  hall  at  position  zbplt  cm. 

From  bottom, 

ACCEL,  IN  Z direction  ( M/SEC**Z) 

Z CUUROS,  OF  rings  of  charge  in  GREEN'S  FCN,  (C“) 

Z POSITION  OF  followed  part,  at  given  time  (CM) 

FUR  10PT(1S)«1,  Z(1  IS  LOwER  Z vALuE  OF  tmE  FINITE  VOLUME 
element  of  charge  in  the  GREEN'S  FCN  SOURCE  TER“ 

HOTTilM  OF  CVL  (CM)  cvl,e*tends  from  0 to  L 

TOP  OF  CVL,  (CM)  (SEE  L) 


Z. POSITION  OF 
Z. POSITION  OF 
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